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Abstract 

We investigate transverse electromagnetic waves propagating in a plasma influenced by the gravitational 

q ■ field of the Schwarzschild-de Sitter black hole. Applying 3+1 spacetime split we derive the relativistic 

O" 1 

two-fluid equations to take account of gravitational effects due to the event horizon and describe the set 

*-i '■ 

of simultaneous linear equations for the perturbations. We use a local approximation to investigate the 

one-dimensional radial propagation of Alfven and high frequency electromagnetic waves. We derive the 

^- ' dispersion relation for these waves and solve it for the wave number k numerically. 
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1 Introduction 

Black holes are still mysterious [T]. Physicists are grappling the theory of black holes, while astronomers 
are searching for real-life examples of black holes in the universe [2J. There still exists no convincing 
observational data which can prove conclusively the existence of black holes in the universe. Black holes, 
however, are not objects of direct observing. They must be observed indirectly through the effects they 
exert on their environment. It is proved that black holes exist on the basis of study of effects which they 
exert on their surroundings. 

Black holes greatly affect the surrounding plasma medium (which is highly magnetized) with their 
enormous gravitational fields. Hence plasma physics in the vicinity of a black hole has become a subject 
of great interest in astrophysics. In the immediate neighborhood of a black hole general relativity applies. 
It is therefore of interest to formulate plasma physics problems in the context of general relativity. 

A covariant formulation of the theory based on the fluid equations of general relativity and Maxwell's 
equations in curved spacetime has so far proved unproductive because of the curvature of four-dimensional 
spacetime in the region surrounding a black hole. Thorne, Price, and Macdonald (TPM) [3l 2J [6] 
developed a method of a 3 + 1 formulation of general relativity in which the electromagnetic equations 
and the plasma physics at least look somewhat similar to the usual formulations in flat spacetime while 
taking accurate account of general relativistic effects such as curvature. In the TPM formulation, work 
connected with black holes has been facilitated by replacing the holes event horizon with a membrane 
endowed with electric charge, electrical conductivity, and finite temperature and entropy. The membrane 
paradigm is mathematically equivalent to the standard, full general relativistic theory of black holes so 
far as physics outside the event horizon is concerned. But the formulation of all physics in this region 
turns out to be very much simpler than it would be using the standard covariant approach of general 
relativity. 

Buzzi, et.al. [7j [8], using the TPM formalism described a general relativistic version of two-fluid 
formulation of plasma physics and investigated the nature of plasma waves (transverse waves in [7] , and 
longitudinal waves in [8]) near the horizon of the Schwarzschild black hole. In this paper we apply 
the formalism of Buzzi et.al. [7] to investigate the transverse electromagnetic waves propagating in a 
plasma close to the Schwarzschild black hole in the de Sitter (dS) space, called the Schwarzschild-de 
Sitter (SdS) black hole. The study is interesting because of the presence of cosmological constant (A). 
There has been a renewed interest in cosmological constant as it is found to be present in the inflationary 
scenario of the early universe. In these scenario the universe undergoes a stage where it is geometrically 
similar to dS space |SJ. Among other things inflation has led to the cold dark matter. According to the 
cold dark matter theory, the bulk of the dark matter is in the form of slowly moving particles (axions 
or neutralinos) . If the cold dark matter theory proves correct, it would shed light on the unification 
of forces [10\ lllj . Comprehensive review works of the cosmological constant or dark energy, including 
the observational evidence of it and the problems associated with it, have been done by many authors 
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In recent years, considerable attention has been concentrated on the study of black holes in dS spaces. 
This is motivated, basically, by the following two aspects: first, several different types of astrophysical 
observations indicate that our universe is in a phase of accelerating expansion [2TJ [22J HH HH H3 HH1 
[23[2S[2S[3Ol[3TJ[3![33[10[3^ Associated with this 

acceleration is a positive cosmological constant. Our universe therefore might approach a dS phase 
in the far future. Second, similar as the AdS/CFT correspondence, an interesting proposal, the so-called 
dS/CFT correspondence, has been suggested that there is a dual relation between quantum gravity on a 
dS space and Euclidean conformal field theory (CFT) on a boundary of dS space [43 HH [49l [50l [51j [52] . 
In view of the above reasons, our study of transverse wave propagation in relativistic two-fluid plasma 
in the environment close to the event horizon of the SdS black hole is interesting. The result we have 
obtained reduces to that of the Schwarzschild black hole [7] when the cosmological constant vanishes. 

This paper is organized as follows. In section [2j we summarize the 3+1 formulation of general rela- 
tivity. In section [3j we describe the nonlinear two-fluid equations expressing continuity and conservation 
of energy and momentum. The two-fluids are coupled together through Maxwell's equations for the elec- 
tromagnetic fields. In section [4[ we restrict one-dimensional wave propagation in the radial z (Rindler 
coordinate system) direction, and linearize the equations for wave propagation in section [5] by giving a 
small perturbation to fields and fluid parameters. We express the derivatives of the unperturbed quan- 
tities with respect to z. In section we discuss the local or mean-field approximation used to obtain 
numerical solutions for the wave dispersion relations. We describe the dispersion relation for the trans- 
verse waves. We give the numerical procedure for determining the roots of the dispersion relation in 
section [7] In section [5] we present the numerical solutions for the wave number k. Finally, in section [5J 
we present our remarks. We use units G = c = ks = 1- 

2 Formalism in 3+1 Spacetime 

In this section we apply TPM formulation to split the Schwarzschild-de Sitter black hole spacetime, which 
is the solution of Einstein equations with a positive A(= 3/£ 2 ) term corresponding to a vacuum state 
spherical symmetric configuration. The metric of the spacetime is asymptotically de Sitter and has the 
form 

ds = g jiV dx^dx v 

= -A 2 dt 2 + -^dr 2 + r 2 (d9 2 + $m 2 0dy 2 ), (1) 



where the metric function is 



.2 2Af r 2 

A 2 = l ^, 2 

r t A 



M being the mass of the black hole and the coordinates are defined such that — oo < t < 00, r > 0, 
< < n, and < < 2-7T. At large r, the metric ([1]) tends to the dS space limit. The explicit dS case 
is obtained by setting M — while the explicit Schwarzchild case is obtained by taking the limit £ —>■ oo. 
When £ 2 is replace by — £ 2 , the metric ([T]) describes an interesting nonrotating AdS black hole called the 



Schwarzschild Anti-de Sitter (SAdS) black hole. The AdS black holes play a major role in the AdS/CFT 
correspondence [5J2 021 ES ES] and they have also received interest in the context of brane- world scenarios 
based on the setup of Randall and Sundrum [STJ EB] ■ 

The horizons of the SdS black hole are located at the real positive roots of A 2 (r) = -pz{r — rh)(r — 
r c )(r_ — r) = 0, and there are more than one horizon if < 3 < 1/27 where 3 = M 2 /£ 2 . The black hole 
(event) horizon ry L and the cosmological horizon r c are located, respectively, at 

2M IT + jj 

rh = -=cos^— , (3) 

2M n-ip 

=cos— — , (4) 



where 

V> = cos _1 (3a/33). (5) 

The negative root 

2M 4, 

cos— (6) 
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has no physical meaning. In the limit 3 — > 0, one finds that fh — > 2M and r c — > £, and it is obvious that 
t c > r h, i-e., the event horizon is the smallest positive root. The spacetime is dynamic for r < Th and 
r > r c . The two horizons coincide: r/, = r c = 3M (extremal), when 3 = 1/27, and the spacetime then 
becomes the well known Nariai spacetime. Expanding r^ in terms of M with S < 1/27, we obtain 

/ AM 2 \ 

n^2A/(^i + — + ...J, (7) 

that is, the event horizon of the SdS black hole is greater than the Schwarzschild event horizon, rs c h = 2M. 
For 3 > 1/27, the spacetime is dynamic for all r > 0, that is, the metric (JT|) then represents not a black 
hole but an unphysical naked singularity at r = 0. 

The hypersurfaces of constant universal time t define an absolute three-dimensional space described 
by the metric 



1 

A 2 



ds 2 = gijdafda? = -r^dr 2 + r 2 (d6 2 + sin 2 ddcp 2 ), (8) 



where the indices i,j range over 1, 2, 3 and refer to coordinates in absolute space. The fiducial observers 
(FIDOs), the observers at rest with respect to this absolute space, measure their proper time r using 
clocks that they carry with them and make local measurements of physical quantities. Then all their 
measured quantities are defined as FIDO locally measured quantities and all rates measured by them are 
measured using FIDO proper time. The FIDOs use a local Cartesian coordinate system with unit basis 
vectors tangent to the coordinate lines 

d_ ^ ld_ ^ _L 8_ 
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For a spacetime viewpoint this set of orthonormal vectors also includes the basis vector for the time 
coordinate, given by 

d ld fi^ 

e ° = J? = ad-f (10) 



where a is the lapse function (or redshift factor) defined by 

°W S I = ( 1 -T-?) ■ (u) 

The FIDO proper time r functions as a local laboratory time, where the FIDOs have the role of "local 
laboratories" . The coordinate time t slices spacetime in the way that the FIDOs would do physically. 

The lapse function a, which governs the ticking rates of clocks and redshifts, plays the role of a 
gravitational potential. It also provides the gravitational acceleration felt by a FIDO [31 HI [5j [6] : 

1 (M r\ 
a=-Vtoa— -(^- ? Je,. (12) 

The rate of change of any scalar physical quantity or any three-dimensional vector or tensor, as measured 
by a FIDO, is defined by the derivative 

D fid 
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v being the fluid velocity measured locally by a FIDO. 

3 Two-Fluid Equations in 3+1 Formalism 

We consider two-component plasma such as an electron-positron plasma or electron-ion plasma. The 
continuity equation for each of the fluid species in the 3+1 notation is 

d 

— {j 8 n s ) + V- (aj 8 n s v s ) = 0, (14) 

where s is 1 for electrons and 2 for positrons (or ions). For a perfect relativistic fluid of species s in 
three-dimensions, the energy density e s , the momentum density S s , and stress-energy tensor W^ k are 
given by 

e s = 7 , 2 (£, + Psvl), S s = 7 2 (e s + P s )v s , Wi k = 7s 2 (e, + Ps)v\v k s + P s g> k , (15) 

where v s is the fluid velocity, n s is the number density, P s is the pressure, and £ s is the total energy 
density defined by 

e s — m s n s + P s /{j g - 1), (16) 

the gas constant j g being 4/3 as T — > oo and 5/3 as T — v 0. 

The ion temperature profile is closely adiabatic and it approaches 10 12 K near the horizon [59]. Far 
from the event horizon, electron (positron) temperatures are essentially equal to the ion temperatures. 
But the electrons closer to the horizon are progressively cooled to about 10 8 — 10 9 K by mechanisms such 
as multiple Compton scattering and synchrotron radiation. Using the conservation of entropy one can 
express the equation of state in the form: 

DIP ^ 0, (17) 



where D/Dt = (l/a)d/dt + v s • V. For a relativistic fluid, the full equation of state as measured in the 
fluid's rest frame, is given by [60| [61]: 

P s iH ( 2 iy (im s n s /P s ) 



e = m s n s + m s n s 



m s n s 



H^ >(im s n s /P s 



(18) 



where the H\ ' (x) are Hankel functions. 

The fluid quantities in (fT5|) take the following form in the electromagnetic field: 

e s = -L(E 2 +B 2 ), S s = iExB, 

87T 47T 

W 3 s k = — (E 2 +B 2 W fe ~— (E J E k + B j B k ). (19) 

87T 4"7T 

Energy and momentum conservation equations are given by [31 SI [5] 

-JU = -V-S s + 2a-S s , (20) 

a at 

-^-S s = Cs a--V-(aW s ). (21) 

a at a 

When the two-fluid plasma couples to the electromagnetic fields, Maxwell's equations take the following 
3+1 form: 

VB = 0, (22) 

V-E = 47R7, (23) 

— = -Vx(aE), (24) 

<9E 

— = V x (aB) - 47raJ, (25) 

where the charge and current densities are respectively defined by 

Q" = y^ 7sqsn s , J = ^2 7 s q s n s v s . (26) 

S S 

Using HHJ) and ([22] US]), the energy and momentum conservation equations ([20]) and (f2T]) can be 
rewritten for each species s in the form 

1 Pi 1 £i 

-^-P 8 - -ut[7 s 2 (£, + Ps)] ~ V • [ 7s 2 (e, + P s )v s ] + 7,3 s n s E ■ v s + 2 7 2 (e s + P s )a ■ v. = 0, (27) 

a at a at 

and 

7s 2 (e s + P.) f -— + v s • V J v s + VP S - j s q s n s (E + v s x B) 

+v s h s q s n s E • v s + -foPs) + 7. 2 (e.s + P«)[v.(v a • a) - a] = 0, (28) 

respectively. Although these equations are valid in a FIDO frame, they reduce for a = 1 to the corre- 
sponding special relativistic equations [62j which are valid in a frame in which both fluids are at rest. 
The transformation from the FIDO frame to the comoving (fluid) frame involves a boost velocity, which 
is simply the freefall velocity onto the black hole, given by 

Uff = (l-a 2 )i (29) 



The relativistic Lorentz factor then becomes 7boost = (1 — W I) -1 ^ 2 = l/ a - 

The two-fluid equations in SdS coordinates form the basis of the numerical procedure for solving the 
linear two-fluid equations; but they cannot be evaluated analytically. However, the Rindler coordinate 
system, in which space is locally Cartesian, provides a good approximation to the SdS metric near 
the event horizon. Without the complication of explicitly curved spatial three-geometries, the essential 
features of the horizon and the 3+1 split are preserved. 

The SdS metric (TT|) in Rindler coordinates is approximated by 

ds 2 = -a 2 dt 2 + dx 2 + dy 2 + dz 2 , (30) 

where 

x = r h ( 9 - — J , y — r h ip, z = 2r h A. (31) 

The standard lapse function in Rindler coordinates becomes a = z/2r^, where r^ is the event horizon of 
the black hole. 

4 Radial Wave Propagation in One-Dimension 

We consider one-dimensional wave propagation in the radial z direction. Introducing the complex vari- 
ables 

v sz (z,t) = u s (z,t), v s (z,t) = v sx (z,t) +iv sy (z,t), 
B(z,t)=B x (z,t)+iB y (z,t), E{z 1 t)=E x {z,t)+iE y {z,t) 1 (32) 

and setting 

v sx B y - v sy B x = -(v s B* ~ V*B), 

v sx E y - v sy E x = t;( v s e * - v *s E )> (33) 

where the * denotes the complex conjugate, one could write the continuity equation (|14[) in the form 

d . . d 
— (7 s n s ) + -^-(aj s n s u s ) = 0, (34) 

and Poisson's equation (|2"3")l as 

BE 

-— ^ = 47r(gini7i + 92^272)- (35) 

Bz 

From the e^ and e^ components of (f24l) and (|25|) . one could obtain 

^f - ->&-)*■ 

BE / B \ 

i— = -a I aj B- i47rea(72n 2 w 2 -71^1^1)- (37) 

Differentiating (|3"T|) with respect to t and using (|36[) . we obtain 

2 B 2 3a 8 B 2 1 \ d . 

E — inea— (^272^2 - riijiVi). (38) 
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The transverse component of the momentum conservation equation is obtained from the e^ and e^ 
components of (|28|) as follows: 



Dv s ( 1 dP s \ 
Psyy 1 = q s n s ^s(E - iv 8 B z + \u s B) - u s v s p s a - v s I q s n s j s E ■ v s -\ -— ^ J , (39) 

where 

E • v s = -(Ev* + E*v s ) + E z u s , 

and p s is the total energy density defined by 

p s = 7s 2 (e s + P„) = ~/ 2 s (m s n s + T g P s ) (40) 

with T g = 7s/(7 s - 1). 

5 Linearization 

We apply perturbation method to linearize the equations derived in the preceding section. We introduce 
the quantities 

u s (z,t) = u os (z) + 5u s (z,t), v s (z,t) = Sv s (z,t), 

n s (z,t) = n os (z) + 6n s (z,t), P s (z,t) = P os (z) + SP s (z,t), 

p s {z,t) = p os {z) + 5p s {z,t) 1 E(z,t) =5E(z,t), 

B z (z,t) = B {z) + SB z (z,t), B(z,t) = 5B(z,t), (41) 

where an applied magnetic field has been chosen to lie along the radial e z direction. The relativistic 
Lorentz factor is also linearized such that 

Is = 7os + Sj s , where j os = (l - u^) 5 , 6f s = ^ s u os ■ 8u s . (42) 

Near the event horizon the unperturbed radial velocity for each species as measured by a FIDO along e z 
is assumed to be the freefall velocity so that 

u os (z) = v s {z) = [l-a 2 {z)}?. (43) 

It follows, from the continuity equation (|3"4")) . that 

r 2 aj os n os u os = const. = r\a h ^ h n h u h , 

where the values with a subscript h are the limiting values at the event horizon. The freefall velocity at 
the event horizon becomes unity so that Uh = 1. Since u os — vs, 7 os = l/cc; and hence aj os — othlh = 1. 
Also, because vs = {i"h/i") 1 Xj the number density for each species can be written as follows: 

n os (z) = — n hs v%{z), (44) 

A 



where 

/ t- 3 \ 2 r /i i m 2 

(45) 



rhr c r- 



, 1 1 

1 + r h — + — 



When i 2 — » oo, (|44|) reduce to the Schwarzschild result |7j. The equation of state ([17]) and (J44j) lead to 
write the unperturbed pressure in terms of the freefall velocity as 



Pos(z) 



37 9 



X 



4 1; 



■PhsV^iz) 



Since P os = kBn os T OSl then with fcg = 1 the temperature profile is 

1 _ „3( 7g -l) 



Tos(z) 



-Th s Va 



(z). 



(46) 



(47) 



The unperturbed magnetic field, chosen to be in the radial z direction, is parallel to the infall fluid 
velocity u os (z)e^ for each fluid. It thus does not experience effects of spatial curvature along with the 
infall fluid velocity. The magnetic field depends on the radial coordinate r because of flux conservation. 
Since V • B = 0, it follows that r 2 B (r) — const., from which one could obtain the unperturbed magnetic 
field in terms of the freefall velocity as 

f 



BJz) 



X 1 



{ B h v%{z). 



Since 



we have 





dvg 

dz 


a 1 
2r h vs ' 




du os 
dz 


a 1 
2r h v s ' 


dB Q 

dz 


4a B 

2r h v% ' 


dn os 


3a n os 


dP os 


3a JgP OS 


A-, 


J*, 


9*.. ,,2 



(48) 
(49) 



(50) 

UZ ,677, V S az Al'h Va 

When the linearized variables from (J4TJ) and (j42|) are substituted into the continuity equation and 
products of perturbation terms are neglected, it follows that 



lo 



d_ 

di 



u os a 



0_ 
dz 



II: 



OS 2 ^^0£ 

2rh ° s dz 



-noslos 



i) 



i) 



1 



dt dz 2rh V n nf! dz 



Sn s - 
1 dn„ 



8 

dz 



+ Sloshes 



1 

2n 

du 



dz 



yl'os^fosV'os ) 

Su s = 0. 



The conservation of entropy, (|17[) . on linearization, gives 

l^os 

Then from the total energy density, ([40]) . it follows that 

Pos f, , llslgPt 



5n s 



Sp s 



1 



Sn s + 2u os -f os p os 5u s 



(51) 



(52) 



(53) 



T^os \ Pos 

where p os — j 2 s (m s n os +T g P os ). Linearizing the transverse part of the momentum conservation equation 
differentiating it with respect to t, and then substituting from (|36|) . we derive 

d d u os iaq s "f 0S n 0S B \ dSv s aq s "/ os n os ( d d u. 



au os 



dz dt 2rt 



' dz dt 2rh pos J dt p os 

When linearized, Poisson's equation ([55)) and ([55]) . respectively, give 

dft F 1 

47re(n o2 7o2 - n ol j i) + ^■Ke(^ o2 8n 2 - JoiSrii) + 47re(n o2 u o2 7o2 (5u 2 - n itt i7oi^i)) 



dz 



d 2 



3a d d 2 



1 



a 



dz 2 2r h dz dt 2 {2r h f 



SE = Airea ( n 27o2 



dSv 2 
~~dT 



noiloi- 



d8v\ 

~dT 



(54) 

(55) 
(56) 



6 Dispersion Relation 

We restrict our consideration to effects on a local scale for which the distance from the horizon does 
not vary significantly. We apply a local (or mean-field) approximation for the lapse function a{z) and 
hence for the equilibrium fields and fluid quantities. If the plasma is situated relatively close to the event 
horizon, a 2 -C 1, then a relatively small change in distance z will make a significant difference to the 
magnitude of a. Thus it is important to choose a sufficiently small range in z so that a(z) does not vary 
much. 

We consider thin layers in the e^ direction, each layer with its own a , where a is some mean value 
of a within a particular layer. By considering a large number of layers within a chosen range of a values, 
a more complete picture can then be built up. However, such a local approximation imposes a restriction 
on the magnitude of the wavelength and, hence, on the wave number k. It is assumed that the wavelength 
is small in comparison with the range over which the equilibrium quantities change significantly. Then 
the wavelength must be smaller in magnitude than the scale of the gradient of the lapse function a, i.e., 

or, equivalcntly, 



A < I — ] = 2r h ~ C5.896 x 10 a cm, 



2-7T 

k> — ~C 1-067 x KT 5 cm-\ 1 < C < 1-5, 
2r h 

for a black hole of mass ~ IMq, The value for £ = 1.5 corresponds to the extremal SdS black hole, while 
the value for £ = 1 corresponds to the Schwarzschild black hole. 

The hydrodynamical approach used in this work has some disadvantages, such as it is essentially 
a bulk, fluid approach. So the microscopic behavior of the two-fluid plasma is treated in a somewhat 
approximate manner via the equation of state. It means that the results are really only strictly valid in the 
long wavelength limit. However, the restriction, imposed by the local approximation, on the wavelength 
is not too severe and permits the consideration of intermediate to long wavelengths so that the small k 
limit is still valid. 

The unperturbed field and fluid quantities are assumed not to be constant with respect to a(z). Then, 
using the local approximation for a, the derivatives of the equilibrium quantities can be evaluated at each 
layer for a given a . In the local approximation for a, a ~ a is valid within a particular layer. Hence, 
the unperturbed fields and fluid quantities and their derivatives, which are functions of a, take on their 
corresponding "mean- field" values for a given ao- Then the coefficients in ([5TJ) . (IS"1)) . and (|5"5")l become 
constants within each layer, evaluated at each fixed mean-field value, a — a . So it is possible to Fourier 
transform the equations with respect to z, using plane- wave- type solutions for the perturbations of the 
form ~ e i ( kz -^ Jt ) f or each a layer. 

When Fourier transformed, (I5"il) and (1551) turn out to be 



uj a ku os - lu + 1 dv s - ia a ku os -u— - — )dh = 0, (57) 

V 2r ft Pos J Pos V 2r h J 

i'±irea Lj{n 27o2Sv2 - n ij ol 6vi) . . 

~ a k(a k - i3/2r ft ) - ^ - l/{2r h f ' ( ' 

10 



The dispersion relation for the transverse electromagnetic wave modes may be put in the form 



K± (k± ± -^-) 



2 , 1 



a 2 \ Upi{u - u ol K ± ) uj 2 p2 {u - u o2 K±) 1 

° ) (u iK T - uj - a o 0J c i) (u o2 K T - lj + a uj c2 ) ( 



2r h J " (2r h ) 

for either the electron-positron or electron-ion plasma, where K± = a k ± i/2r/j, uj cs — ej 0S n 0S B / p 



and ujp S = \/ ^ e 2 J os n os / Pos- Like the plasma frequency u> ps , the cyclotron frequency uj cs is frame 
independent. Although the fluid quantities are measured in the fluid frame, the field B is measured in 
the FIDO frame. So, the factors of 7 os do not cancel out explicitly. The transformation B — > j sB 
boosts the fluid frame for either fluid and thereby cancels the j os factors. The "+"and "—"denote the 
left L and right R modes, respectively. The dispersion relation for the L mode is obtained by taking 
the complex conjugate of the dispersion relation for the R mode. Both of the two modes have the same 
dispersion relation in the special relativistic case. 

7 Numerical Solution of Modes 

We have to solve the dispersion relation (I59j) in order to determine all the physically meaningful modes 
for the transverse waves. But the dispersion relation is complicated enough, even in the simplest cases for 
the electron-positron plasma where both species are assumed to have the same equilibrium parameters, 
and an analytical solution is cumbersome and unprofitable. We therefore solve the dispersion relation 
numerically and for this purpose we put it in the form of a matrix equation as follows: 

(A - H)X = 0, (60) 

where the eigenvalue is chosen to be the wave number k, the eigenvector X is given by the relevant set 
of perturbations, and I is the identity matrix. 

To write the perturbation equations in an appropriate form, we introduce the following set of dimcn- 
sionless variables: 

w r kc 1 

oj = , k = — , k 



-.„ Su s . Sv s Sn s 

ou s = , v s = , on s = , 

^os t^os Tlos 

For an electron-positron plasma, aj pl = w p2 and w c i = lu c2 , since the choice of input parameters is the 

same for each fluid. Hence, a;* is defined as 

lu c Alfven modes, 

(62) 
(2cjp + oj 2 ) 2 high frequency modes, 

where io v = Juj v \lo v2 and uj c — ^/uj c iuj c2 . However, for the case of an electron-ion plasma, the plasma 

frequency and the cyclotron frequency are different for each fluid, and so the choice of w* is a more 

complicated matter. For simplicity, it is assumed that 

77=(u>ci + w c2) 5 Alfven modes, 
w* = ^ (63) 

(ui^i + oj 2 2 )i high frequency modes, 

11 



where uj^ s = {2loI s +uj^ s ). These chosen values for u>* reduce to the special relativistic cutoffs in the zero 
gravity limit for an electron-positron plasma. The cutoffs, in the special relativistic case, are determined 
by the dispersion relation in view of the fact that the solutions to the dispersion relation are physical 
(i.e., Re(fc) > 0) only for certain frequency regimes. Since the dispersion relations cannot be handled 
analytically, it is difficult to determine the cutoffs in the case including gravity. Other similar combinations 
for w* make only a difference of a scale factor to the form of the results as a;* is really. 
We write the dimensionless eigenvector for the transverse set of equations as 



X 



transverse 



Svi 

Sv 2 

SB 
5E 

When linearized and Fourier transformed, equations (I36|) and (I37|) turn out to be 



(64) 



k —)SE+—SB = 0, (65) 



and 

i^ 



5E = [k — SB + 47re(7o2n-o2<5w2 - "/oinoiSvi). (66) 

a \ 2r h a J 

Using ([61]) . we write ([57]) . (I65| . and (|66|) in the dimensionless form as follows: 



kSv* — 



I _ /&\ _^_ _ jM Bs + /«.\ J^ s§ j («.\ J^_ sM (67) 

-j S V e / u os lu* a J \ e / u os uj* \ e J u os lo^ 

ikh 
kSE = -iCjSB + —5E, (68) 

a 

kSB = u i-^—dvi-Uo2-^-Sv2 + —SB + iQ6E. (69) 

These are the equations in the required form to be used as input to (|60[) . 

8 Results 

We have carried out the numerical analysis by using the well known MATLAB. We have considered both 
the electron-positron plasma and the electron-ion plasma. The limiting horizon values for the magnetic 
field and the fluid parameters are chosen as follows. For the electron-positron plasma, the horizon values 
are taken to be 

n hs = 10 18 cm~ 3 , T hs = 10 10 K, B h = 3 x 10 6 G, and lg = |. (70) 

For the electron-ion plasma, the ions are essentially nonrelativistic, and the limiting horizon values 
are chosen to be 

n hl = f0 18 cm- 3 , T hl = 10 10 K, n h2 = 10 15 c m - 3 , T h2 = 10 12 K. (71) 

The limiting horizon value for the equilibrium magnetic field has the same value as it has for the electron- 
positron case. The limiting horizon temperature for each species has been taken as derived by Colpi et 
al. [59j from studies of two temperature models of spherical accretion onto black holes. For the limiting 
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Figure 1: Top: Real part of the Alfven growth and damping modes for the electron-positron plasma. 
Center: Imaginary part of the growth mode. Bottom: Imaginary part of the damped mode. 



horizon densities and the limiting horizon field, simply arbitrarily those values have been chosen which 
are consistent with the current ideas. The gas constant and the mass of the black hole have been taken 



MS 



lg 



M = 5M, 



0- 



(72) 



8.1 Alfven Modes 

8.1.1 Electron-Positron Plasma 

For the ultrarelativistic electron-positron plasma in the special relativistic case, only one real Alfven 
mode is found to exist 62J, while for the Schwarzschild case there are two Alfven modes [7J. Our study 
in the SdS spacetime shows four modes to exist for electron-positron plasma. The first two modes, shown 
in the Fig. 1 are drawn from the dispersion relation of (f67|) . (|68|) and (|69|) and are complex conjugate 
of each other. They are equivalent to that of Buzzi et al. [7J in this case. The third mode is damped 
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Figure 2: Left: Real part of the Alfven growth mode for the electron-ion plasma. Right: Imaginary part 
of the growth mode. 




Figure 3: Left: Real part of Alfven damped mode for the electron-ion plasma. Right: Imaginary part of 
the damped mode. 

mode similar to the damped mode of Fig. 1 with larger damping rate and the fourth mode is equivalent 
to the third mode with opposite in sense and shows growth. These four modes coalesce into a single 
mode in the special relativistic limit, giving the result of ref. [62j . Since we are using the convention 



ikz 



,i[Re(fc)+iim(fc)]z ^ £ ne o; anl p m g corresponds to Im(fc) > and growth to Im(fc) < 0. 



8.1.2 Electron-ion Plasma 

In this case there exist four modes, two of which are respectively growth and damped as shown in Figs. 
2 and 3. The other two modes, complex conjugate to each other, are equivalent to the electron-positron 
modes shown in Fig. 1 but have larger amount of damping and growth rates. For the first two modes 
the differences in the magnitudes of the uj c \ and lo C 2 apparently lead to take the frequencies from their 
negative (and therefore unphysical) values for the electron-positron case to positive physical values for 
the electron-ion case. These changes are thus because of the difference in mass and density factors as 
between the positrons and ions. It is evident that the growth and damping rates are independent of the 
frequency, but depended only on the distance from the black hole horizon through a . 
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Figure 4: Left: Imaginary part of high frequency damped mode. Right: Imaginary part of high frequency 
mode showing both damping and growth. 

8.2 High Frequency Transverse Modes 
8.2.1 Electron-Positron Plasma 

In this case there are four high frequency electromagnetic modes. Two of these modes, shown in Fig. 4, 
have real parts equivalent to the real part of Fig.l. The left (imaginary) mode in Fig. 4 is damped, while 
the right one is damped for higher frequency and a Q > 0.2 and growth for a < 0.2 and lower frequencies. 
It then appears, at a distance from the horizon corresponding to a < 0.2, that energy is no longer fed 
into wave mode by the gravitational field but begins to be drained from the waves. The third mode, 
shown in Fig. 5, is damped for most of the frequency domain and a > 0.2, but shows growth for lower 
frequencies and a < 0.2. Near to the event horizon, below about a ~ 0.2, it becomes a growth mode for 
all frequencies. The fourth mode, shown in Fig. 6, is similar to third mode but shows a growth for all the 
frequency domain. These four modes reduce to the three modes of Buzzi et al. [7] in the Schwarzschild 
case, while they coalesce into only one purely high frequency mode in the special relativistic case |62j for 
the ultrarelativistic electron-positron plasma. 




Figure 5: Left: Real part of high frequency mode for the electron-positron plasma. Right: Imaginary part 
of high frequency damping and growth mode. 
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Figure 6: Left: Real part of high frequency mode for the electron-positron plasma. Right: Imaginary part 
of high frequency growth mode. 

8.2.2 Electron-Ion Plasma 

Like the electron-positron plasma, the electron-ion plasma admits four high frequency modes. The first 
two modes are similar to the modes for electron-positron case shown in Fig. 4, with a large amount of 
growth and damping rates. The third and fourth modes are similar to the modes shown in Fig. 5 and 
Fig. 6, respectively, with larger amount of damping and growth. The third mode is growth for higher 
frequency and a > 0.1, while it is damped for lower frequency and a — -> 0. The fourth mode is a 
growth mode and is stable for all frequencies and at all distances from the event horizon. Unlike the 
corresponding Alfven modes, the growth and decay rates obviously depend on frequency. The four modes 
reduce to three modes in the Schwarzschild case [7] . 

9 Concluding Remarks 

The study of this paper concerns exclusively the investigation, within the local approximation, of Alfven 
and high frequency transverse electromagnetic waves in a two-plasma surrounding the Schwarzschild black 
hole in the de Sitter space. Using a local approximation, we have derived the dispersion relations for the 
Alfven and high frequency electromagnetic waves, and solved it numerically as an analytical solution is 
unprofitable and not illuminating. In the limit of zero gravity our results reduce to the special relativistic 
results |62j where only one purely real mode was found to exist for both the Alfven and high frequency 
electromagnetic waves. In contrast to the special relativistic case, new modes (damped or growth) arise 
in the black hole spacetime. This is because of the black hole's gravitational field. 

The damping and growth rates for the electron-positron plasma are smaller in general, by several 
orders of magnitude, compared with the real components of the wave number. However, there exist 
modes for the electron-ion plasma for which the damping and growth rates are significant for the Alfven 
waves, in particular. For the Alfven waves, the damping and growth rates are obviously frequency 
independent, but are dependent on the radial distance from the horizon as denoted by the mean value of 
the lapse function a . This is of course not the case for the high frequency waves for which the rate of 
damping or growth is dependent on both frequency and radial distance from the horizon. 
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Presence of damped modes demonstrates, at least in this approximation, that energy is being drained 
from some of the waves by the gravitational field. Since the majority of the modes are growth rates, it 
indicates that the gravitational field is, in fact, feeding energy into the waves. 

In the limit i 2 — > oo our study provides the results for the Schwarzschild black hole case [7] , while for 
M = the results reduce to those of the pure de Sitter space [63]. When M 2 = 27£ 2 , the results go for the 
well known Nariai spacetime. If £ 2 is replaced by — £ 2 , our study gives results for the interesting AdS black 
holes which have received interest in the AdS/CFT correspondence [53, 54, 55, [5B] as well as in the context 
of brane- world scenarios [57, 58]. In fact, development in string/M theory have greatly stimulated the 
study of black holes in AdS space. Moreover, recent astrophysical observations indicate that our universe 
is in a phase of accelerating expansion associated with which is a positive cosmological constant, and 
the universe therefore might approach a de Sitter phase in the far future [21] [24J [25[ [BH [65[ [66, 67J. In 
view of these reasons, aspects of black holes in the de Sitter space might be of interest, and our study 
of two-fluid plasma in this paper near the event horizon of the Schwarzschild-de Sitter black hole is thus 
well motivated. 
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